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ABSTRACT 

We study a sample of 883 sources detected in a deep Very Large Array survey at 1.4 GHz in 
the Extended Chandra Deep Field South. The paper focuses on the identification of their optical 
and infrared (IR) counterparts. We use a likelihood ratio technique that is particularly useful when 
dealing with deep optical images to minimize the number of spurious associations. We find a reliable 
counterpart for 95% of our radio sources. Most of the counterparts (74%) are detected at optical 
wavelengths, but there is a significant fraction (21%) only detectable in the IR. Combining newly 
acquired optical spectra with data from the literature we are able to assign a redshift to 81% of the 
identified radio sources (37% spectroscopic). We also investigate the X-ray properties of the radio 
sources using the Chandra 4 Ms and 250 ks observations. In particular, we use a stacking technique 
to derive the average properties of radio objects undetected in the Chandra images. The results of 
our analysis are collected in a new catalog containing the position of the optical/IR counterpart, the 
redshift information and the X-ray fluxes. It is the deepest multi- wavelength catalog of radio sources, 
which will be used for future study of this galaxy population. 

Subject headings: cosmology: observations - galaxies: active galaxies: starburst - radio continuum: 
galaxies 



1. INTRODUCTION 

Deep radio observations provide a powerful opportu- 
nity to investigate the high redshift Universe. More- 
over, since radio observations are almost unaffected by 
dust extinction, they allow us to observe objects, which 
are heavily obscured in other bands. While bright radio 
sources are mostly powerful radio galaxies and radio-loud 
(RL) active galactic nuclei (AGN), at lower flux densities 
we observe an increasing fraction of star- forming galaxies 
(SFG) and radio-quiet (RQ) AGN (e.g., iPadovani et al.l 
2009, and references therein). 

Source classification of deep radio surveys is not easy 
and requires multi-frequency data. This approach was 
adopted in a series of papers, which studied a radio se- 
lected sample o f 266 objects in th e Chandra Deep Field 
South (CDFS; [K ellerma nn et al.l [200 8: Mai nieri et aD 
120081 : iTozzi et al.l 120091 : IPadovani etaTI I2009L I2011D . 
Combining the information from different wavelengths, 
these authors were able first to classify the sources as 
SFG, RQ AGN and RL AGN and then to study the prop- 
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erties and the evolution of the different classes separately. 
The promising results of this work have encouraged us 
to apply it to a new radio catalog (N. Miller et al. 2012, 
in preparation) that reaches a lower flux density limit 
and has a more uniform coverage of the Extended- CDFS 
(E-CDFS). As a consequence, we have three times more 
objects, with most of the new sources in the sub-mJy 
regime (at 1.4 GHz). 

This paper focuses on the identification of the optical 
and IR counterparts of the radio sources. Our main goal 
is to assign a redshift to the radio sources and to associate 
them with the correct photometry. This information will 
be then used in future papers to classify the sources and 
study their evolutionary properties. As mentioned above, 
a faint radio selected sample includes sources of widely 
different nature: SFG with a blue stellar population to- 
gether with radio galaxies commonly hosted in redder 
objects, and obscured AGN together with bright unob- 
scured quasars. For this reason, it is important to con- 
sider a large wavelength range, from the ultraviolet to 
the mid-infrared (MIR). Faint radio sources often cor- 
respond to faint optical counterparts. Therefore, deep 
optical observations are needed. This has an impact on 
the methodology that should be adopted in the identifi- 
cation process. 

The structure of the paper is as follows: in Section[2]we 
present the datasets, while in Section [3] we describe the 
likelihood ratio method we used to identify the counter- 
parts of our sources and the results of the identification 
process, including an estimate of the spurious associa- 
tion fraction and a comparison with the cross-correlation 
method. Section [4] discusses the redshift distribution 
(spectroscopic and photometric) of our sample. In Sec- 
tion [5] we deal with the X-ray counterparts of the radio 
sources, while the description of the released catalog is 
given in Section [6] In Section [71 we discuss our results 
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Figure 1. Flux density distribution of the 5cr E-CDFS ra- 
dio i_cata]£g_JsoJid^nej_com to the sample described 
in iKellermani Tet alJ (|2008l ) (dashed line). The 5cr sample 
is about three times larger and the majority of the sources 
have sub-mJy flux densities. 

and report our conclusions. Finally, in Appendix [A] we 
present new redshifts and spectra for the optical coun- 
terparts of 13 Very Large Array (VLA) sources and in 
Appendix iBl we report on some peculiar sources. In this 
paper we use magnitudes in the AB system, if not oth- 
erwise stated, and we assume a cosmology with Ho = 70 
km s _1 Mpc -1 , Q M = 0.3 and ^ A = 0.7. 

2. DATA 

2.1. The radio catalog 

The E-CDFS was observed at 1.4 GHz with the VLA 
between 2007 June and September (Mille ret al.l 12008). 
The mosaic image covers an area of about 34 x 34 arcmin 
with near-uniform sensitivity. The typical rms is 7.4 /iJy 
for a 2.8" x 1.6" beam. The second data release (N. 
Miller et al. 2012, in preparation) provides a new source 
catalog with a 5<r point-source detection limit, for a total 
of 883 sources. We assigned a progressive identification 
number (RID) to the sources ordered by increasing right 
ascension. The flux density distribution of the sample is 
shown in Fig. [TJ where we use the peak flux density or the 
integrated flux density according to the specifications of 
N. Miller et al. (2012, in preparation). The median value 
of the distribution is 58.5 /iJy and the median signal to 
noise ratio (S/N) is 7.6. We note that - 90% of the 
sample has flux density below 1 mJy, a regime where RQ 
AGN and SFG ■ become th e dominant populations (e.g., 
iPadovani et al.ll2QQ9L l2Qllh . A classification of the radio 
sources will be presented in M. Bonzini et al. (2012, in 
preparation). 

2.2. Auxiliary data 

The E-CDFS is one of the most studied patches of the 
sky and has been observed in many wavebands. As we 
will discuss in the following sections, this wealth of data is 
crucial to select the correct counterpart of a radio source. 
Here we describe the large amount of optical and IR data 
used in this work. We considered a total of ten catalogs. 
The complete list is reported in Table [T] together with 



some basic information: the instrument used (column 
2), the effective wavelength (column 3), the typical point 
spread function (PSF) (column 4), the 5<r AB magnitude 
limit (column 5), and the total area covered (column 6). 
The latter is also shown in Fig. [2j where the footprint 
of each mosaic image is plotted over the VLA image. 
For details on the different data sets, we refer to the pa- 
pers listed in column 7. We divide the auxiliary catalogs 
in three groups according to their selection band: op- 
tical, near-infrared (NIR) and mid-infrared (MIR). The 
first group includes U-VIMOS, v-GEMS, R-WFI, and z- 
GOODS. We note that the Wide Field Imager (WFI) ob- 
servations are the only optical images covering the whole 
VLA area. Therefore, even if they are shallower than 
the others and with a lower spatial resolution, they were 
crucial in the identification process. The U-VIMOS cat- 
alog has been produced by us using the SExtractor soft- 
ware (|Bertin fc Arnou ts 1996) from the original images. 
In the NIR, we used the H-GNS, H-SOFI, Ks-MUSYC 
and Ks-ISAAC catalogs. The H-GNS data, consisting of 
30 pointed observations, cover a very small area of the 
E-CDFS but have a better resolution compared to the 
ground-based observations. At longer wavelengths, the 
E-CDFS was mapped with the Spitzer Space Telescope 
as part of the SIMPLE and FIDEL surveys. These data 
are particularly useful to identify high redshift sources 
(see Sec. 14.41 for details). 

3. COUNTERPART IDENTIFICATION METHOD 

3.1. Likelihood ratio technique 

The first step in the identification process consisted 
of registering each auxiliary catalog to the astrometric 
frame of the radio image, correcting for the median off- 
sets between the radio and the auxiliary catalogs. An 
average number of 400 sources was used to perform this 
registration and the typical median offset is 0.2". As 
already mentioned, a simple cross-correlation method, 
where the counterpart is selected as the closest object to 
the radio source given a threshold matching radius, can 
lead to a large number of spurious association when deal- 
ing with deep optical images. Therefore, we adopted a 
likelihood ratio techniqu e (e.g.. iSutherland fc Saunders I 
1992: ICiliegi et aLll2003[ ). This method allows us to take 
into account not only the position of the counterpart, 
but also the background source magnitude distribution 
and the presence of multiple possible counterparts for the 
same radio source. Here we briefly descri be this tech- 
nique following the formalism described in ICiliegi et al.l 
(2003). It consists of three main steps: 

(a) Compute the surface density of background sources 
n(m) as a function of magnitude m. 

(b) Evaluate the likelihood ratio (LR) for each possible 
counterpart. 

(c) Compute the reliability (rel) of each association. 

(a) The magnitude distribution of background sources 
is obtained by counting all the sources within a 30" 
radius around each radio object and dividing them in 
magnitude bins (Am = 0.5). The size of the search- 
ing radius is set to contain, on average, just one radio 
source and a substantial number (>100) of background 
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Table 1 

Auxiliary catalogs used for the identification of the radio sources counterparts. A description of the columns is given in Section [221 





(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 




Catalog 


Instrument 


X eff 


bb b WnM 


AB mag 


Area 


References 








(/mi) 


(arcsec) 


[ocr limit ) 


(arcmin ) 






U-VIMOS 


VIMOS/ESO VLT 


0.390 


0.2 


28.0 


~800 


Nonino et al. (2009) 


Optical 


v-GEMS 


ACS/ HST 


0.578 


0.1 


28.5 


~800 


Rix et al. (2004) 
















Caldwell et al. (2008) 




R-WFI 


WFI/ESO 2.2m 


0.654 


0.8 


25.5 


>1100 


Giavalisco et al. (2004) 




z-GOODS 


ACS/HST 


0.912 


0.1 


28.2 


~160 


Dickinson et al. (2003) 
















Giavalisco et al. (2004f Name 




H-GNS 


NICMOS/HST 


1.607 


0.2 


26.5 


-43 


Conselice et al. (2011) 


NIR 


H-SOFI 


SOFI/ESO NTT 


1.636 


0.55 


22.0 


~800 


Olsen et al. (2006) 


Ks-ISAAC 


ISAAC/ESO VLT 


2.745 


0.4 


24.7 


~131 


Retzlaff et al. (2010) 




Ks-MUSYC 


ISPI/CTIO 4m 


3.323 


0.3 


22.3 


~900 


Tavlor et al. (2009) 


MIR 


IRAC-SIMPLE 


IRAC/Spitzer 


3.507-4.436 c 


1.7 


23.8-23.6 


>1100 


Damen et al. (2011) 




24um-FIDEL 


MIPS/ Spitzer 


23.209 


6 


20.2 


>1100 


Dickinson & FIDEL team (2007) 



used in the text to refer to a specific catalog. 
b Filter effective wavelength. 

c The catalog is obtained from a combined image of 3.6/im and 4.5/im IRAC bands. 
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Figure 2. Multi- wavelength coverages overplotted on the VLA image. Left: Optical catalogs (from outside): v-GEMS (solid 
line), U-VIMOS (dot-dashed line), and z-GOODS (dashed line). The R-WFI catalog coverage exceeds the VLA image and 
therefore is not plotted. Middle: NIR catalogs (from outside): Ks-MUSYC (dot-dashed line), H-SOFI (solid line), Ks-ISAAC 
(dashed line) and H-GNS (small squares). Right: MIR catalogs: IRAC-SIMPLE (dashed line) and 24um-FIDEL (solid line). 
The dot-dashed line encloses the area with photometric redshift catalogs coverage (see Section 

sources for the deep optical catalogs. The surface den- 
sity of background objects n(m) is then computed divid- 
ing the obtained distribution by the total searching area 

(jT X (30 ) X N rac [io sources)- 

(b) The likelihood ratio for a counterpart candidate 
is denned as the ratio between the probability that the 
source is the correct identification and the corresponding 
probability for an unrelated background source. There- 
fore we compute LR as: 



(Jo. 



eri 



LR 



q(m)f(r) 
n(m) 



(i) 



where f(r) is the probability distribution function of the 
positional errors and q(m) is the expected distribution of 
the counterparts as a function of m. As for f{r) we adopt 
a two-dimensional Gaussian distribution of the form: 



/(r) = J_ exp( _^) 



(2) 



, + er 2 and 

er'j . er a and er$ are the radio positional 
errors given by the beam size (1.6 x 2.8 arcsec) divided 
by two times the S/N ratio of the considered source. To 
account for further uncertainties on the VLA position, 
we added in quadrature 0.1" to the radio positional er- 
ror (N. Miller et al., 2012, in preparation). The aver- 
age positional error er aux for the optical catalog is 0.1" 
and 0.3" for the others, with the exception of H-GNS 
(er aux = 0.2") and 24/im-FIDEL (er aux = 0.6") (see ref- 
erences given in Tabled]). To derive an estimate for q(m), 
we first counted all the objects in the auxiliary catalog 
within a radius of 2" around each radio source. Then, we 
subtracted the distribution of background objects com- 
puted on the same area (n(m) xttx (2") 2 x N ra di sources)- 
The latter is shown in Fig. [3] (dashed line), from left 
to right, for an optical, NIR and MIR catalog, respec- 
tively. The background subtracted distribution, reaZ(ra), 
is plotted in the same figure as a solid line. Finally, we 
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Figure 3. Background source distribution, n(m), multiplied for the searching area (tt x (2") 2 x N ra dio sources) (dashed line) 
and the background subtracted distribution of counterparts, real(m), (solid line) for the v-GEMS catalog (on the left), the 
Ks-MUSYC catalog (in the middle) and the IRAC-SIMPLE catalog (on the right). The dotted dashed line is obtained by 
smoothing the real(m) distribution and it is used to compute q(m). 

normalized the distribution function as: 



q(m) 



real(m) 
J2i real(m), 



x Q 



(3) 



where the sum runs over the total number of objects 
in the real(m) and Q represents the probability that 
the real counterpart is abov e the catalo g; dete ction 
limit. As already v erified by iCiliegi et al.l ([2003) and 
M ainieri et al.l ([2008), the number of identifications and 
the associated reliabilities have a mild dependence on 
Q. Therefore, we adopted a fixed value Q = 0.8 for all 
the auxiliary catalogs as it corresponds to the average 
expected fraction of identifications. Finally, combining 
q(m), f(r) and n(m) according to eq. HJ we computed 
the LR for each source in the auxiliary catalogs. 

(c) The LR does not contain information about the 
possible presence of many counterpart candidates in the 
surrounding of a specific radio source. It is therefore 
useful to define the reliability of each association as: 



relj = 



(LR), 



EiiLR)i + (l-Q) 



(4) 



where the sum is over all the candidate counterparts for 
the same radio source ([Sutherland fc Saunders lll992f ). 

3.2. Identification results 

Following the method described in the previous sec- 
tion, we built a list of possible counterparts for each 
auxiliary catalog. Initially, we set a very low likelihood 
threshold (10 -6 ) to be sure not to lose any counterpart. 
After a careful analysis, we decided to consider as reli- 
able only counterparts with reliability greater than 0.6. 
This threshold ensures that the expected number of spu- 
rious associations is below 5% for each auxiliary cata- 
log (see Section [33]) . and at the same time maximizes 
the number of identified sources. The identification rate 
for each auxiliary catalog is reported in column (3) of 
Table El The number of identified sources is weighted 
by the number of radio sources inside the area covered 
by each survey reported in column (2). We note that 
the number of identifications increases with wavelength, 
from ~ 65% in the optical catalogs up to 87% in the 
MIR. That means that most of the radio sources have a 
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Figure 4. Spitzer /IRAC images of objects RID 018 (top) 
and 866 (bottom) . They represent examples of IRAC sources 
not present in the SIMPLE catalog whose position has been 
manually extracted (see Section \3. 41 for details): RID 018 has 
been observed only in two of the four channels, and RID 866 is 
blended by a nearby source. The position of the radio source 
is marked by a cross. Each image is 10 arcsec on a side. 

counterpart candidate in more than one auxiliary cata- 
log, and that there is a fraction of sources that are not 
detected in the optical but only in the IR. In more detail, 
there are 652 radio sources (74%) that have a counter- 
part in at least one of the four optical catalogs, 76 (9%) 
that have no counterpart in any of the optical catalogs, 
but that are identified in at least one of the NIR catalogs, 
and 111 (12%) that have a counterpart only in the MIR. 
We will refer to these three groups as optical, NIR and 
MIR selected counterparts, respectively. We anticipate 
that they have different redshift distributions, with NIR 
and MIR selected sou rces having on average higher red- 
shift (see Section |4~4|) . High redshift objects, thanks to 
their positive K-correction, are more easily observed in 
the IR than in the optical and this explains the higher 
identification rate observed in the MIR catalogs. 

In Table [3j we report the complete list of the coun- 
terparts of the radio sources (see Section EJ). The coun- 
terpart position is taken from an optical catalog, when 
available, since these observations have the highest spa- 
tial resolution. In particular, we chose the catalog in 
which the counterpart has the highest reliability. Ac- 
cording to this criterion, we selected 104 counterparts 
from U-VIMOS, 150 from v-GEMS, 301 from R-WFI, 
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Table 2 

Counterparts identified in each catalog and spurious association estimate for both the likelihood and cross-correlation 
methods. In column 2 we report the number of radio sources inside the area of the survey (see Fig. [2]). 





LR method cross-correlation 


(1) ^ (2) 
Catalog Radio sources 


(3) (4) (5) 
% counterparts (#) 90% inside a % spurious 
(arcsec) 


(6) (7) 
% counterparts (#) % spurious 


U-VIMOS 540 

Optical V - GEM S 646 
R-WFI 877 
z-GOODS 164 


64% (347) 0.7 5% 
67% (432) 0.7 5% 
68% (600) 0.7 3% 
65% (107) 0.5 4% 


64% (346) 12% 
67% (435) 11% 
67% (587) 6% 
58% (95) 6% 


H-GNS 34 
NIR H-SOFI 523 
Ks-ISAAC 135 
Ks-MUSYC 724 


70% (24) 0.5 <1% 
69% (363) 0.7 2% 
81% (109) 0.7 4% 
77% (556) 0.7 3% 


62% (21) 6% 
65% (339) 1% 
76% (102) 3% 
71% (515) 1% 


MIR IRAC-SIMPLE 858 
24um-FIDEL 878 


87% (746) 0.7 3% 
85% (745) 1.1 4% 


79% (674) 3% 
79% (692) 2% 



Radius within 90% of the counterparts are included. 



Table 3 

Radio information and identification process results. A description of the table content is given in Section [6] (This Table is 
available in its entirety in a machine-readable form in the online journal. A portion in shown here for guidance regarding its 

form and content.) 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) (9) 


(10) 


RID 


RA radio 


Dec radio 


Sr 


S/N 


RA counterpart 


Dec counterpart 


Reliability Distance 


Counterpart 




(J2000) 


(J2000) 


OJy) 




(J2000) 


(J200) 


(arcsec) 


catalog 



360 


3:32:13.09 


-27:43:50.9 


1380.0±26.6 


115.8 










unidentified 


361 


3:32:13.24 


-27:40:43.7 


58.9±6.4 


9.0 


03:32:13.24 


-27:40:43.39 


0.97 


0.2 


R-WFI 


362 


3:32:13.25 


-27:42:41.3 


86.4±6.4 


13.5 


03:32:13.25 


-27:42:40.86 


1.00 


0.2 


z-GOODS 


363 


3:32:13.36 


-27:39:35.2 


48.5±6.9 


7.0 


03:32:13.32 


-27:39:35.03 


1.00 


0.2 


24um-FIDEL 


364 


3:32:13.41 


-27:33:04.9 


48.4±7.8 


6.2 


03:32:13.39 


-27:33:04.93 


0.97 


0.3 


R-WFI 


365 


3:32:13.50 


-27:49:53.1 


44.0±6.4 


6.9 


03:32:13.48 


-27:49:52.82 


1.00 


0.1 


z-GOODS 


366 


3:32:13.61 


-27:34:04.3 


102.6±18.5 


8.3 


03:32:13.58 


-27:34:04.37 


1.00 


0.4 


R-WFI 


367 


3:32:13.65 


-28:01:01.2 


35.7±6.9 


5.1 


03:32:13.62 


-28:01:01.06 


1.00 


0.3 


v-GEMS 


368 


3:32:13.85 


-27:56:00.3 


44.5±6.4 


6.9 


03:32:13.85 


-27:55:59.95 


1.00 


0.2 


24um-FIDEL 


369 


3:32:14.17 


-27:49:10.6 


95.4±6.4 


14.9 


03:32:14.14 


-27:49:10.09 


0.98 


0.4 


U-VIMOS 


370 


3:32:14.46 


-27:45:40.8 


34.9±6.2 


5.6 


03:32:14.43 


-27:45:40.72 


1.00 


0.3 


z-GOODS 


371 


3:32:14.60 


-27:43:05.8 


32.1±6.4 


5.0 


03:32:14.60 


-27:43:06.10 


9.00 


0.3 


manual 


372 


3:32:14.69 


-28:02:20.2 


44.9±7.1 


6.0 


03:32:14.65 


-28:02:19.97 


1.00 


0.4 


v-GEMS 


373 


3:32:14.85 


-27:56:40.9 


109.7±6.5 


16.9 


03:32:14.83 


-27:56:40.49 


1.00 


0.2 


v-GEMS 


374 


3:32:15.17 


-28:05:22.7 


50.8±7.9 


6.3 


03:32:15.14 


-28:05:22.24 


1.00 


0.3 


R-WFI 


375 


3:32:15.34 


-27:50:37.6 


43.1±6.4 


6.7 


03:32:15.32 


-27:50:37.25 


1.00 


0.1 


H-GNS 



and 96 from z-GOODS. If there was no optical coun- 
terpart above the reliability threshold, we used the co- 
ordinates of the most reliable counterpart found in the 
NIR catalogs. This happened for 4 sources from H-GNS, 
24 from H-SOFI, 4 from Ks-ISAAC, and 47 from Ks- 
MUSYC. For the remaining counterparts, we used the 
position from the IRAC-SIMPLE catalog (25 sources) 
and from the 24um-FIDEL one (74 sources). Finally, 
there are 10 radio sources (RID: 18, 19, 36, 371, 430, 
457, 463, 698, 795, and 866) whose counterpart is clearly 
visible in the IRAC images but is not listed in the SIM- 
PLE catalog (or in any other catalog). The reason is 
that, since the SIMPLE catalog is extracted from the 
combined 3.6 and 4.5 fim images, either the source was 
observed only in the first or the second IRAC channel 
and therefore not included in the catalog, or it was not 
deblended from a nearby object (see Fig. 2J. In these 
cases we have extracted the position of the counterpart 



from the IRAC image. 

As a further check, we extracted 10x10 arcsec cutouts 
centered at the radio source position of the images in 
the various bands, to visually inspect the counterpart 
associations. Examples are presented in Fig. [5] where 
the position of the radio source and of its counterpart are 
marked by a cross and a square, respectively. In the left 
panels, radio contours are plotted over a 20x20 arcsec R 
band image, after the latter has been registered on the 
astrometric frame of the radio image. This larger size is 
chosen for a better view of the radio contours. In most 
cases the selected counterpart is clearly visible in one or 
more cutouts. In 12 cases, we found a more convincing 
counterpart and therefore we revisited the associ ation ; 
these cases are discussed in more details in Section [3T4l 

A total of 44 sources are unidentified: most of them are 
either very faint radio sources or lie at the edge of the 
field, where the multi-wavelength coverage is less rich. 
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They are blank fields in all the available images (see, 
e.g., RID 360 in Fig. [5j). We expect only few of them 
to be spurious radio detections since the radio catalog is 
based on a mosaic image and therefore each object was 
observed by more than one pointing (N. Miller et al., 
2012, in preparation). If an object were spurious and 
due to instrumental effects (e.g., a sidelobe of a nearby 
bright source) it would not be the same in each point- 
ing. Similarly, if it were just a chance noise spike you 
would expect to see it in only one pointing. Therefore 
we perform source fits in the individual pointings for the 
unidentified radio sources with low S/N and we believe 
that the radio detections are real, with one or two possi- 
ble exceptions. 

In summary, we found a reliable counterpart for 839 
out of 883 radio sources (95%). 

3.3. Multiple component radio systems 

A multi-wavelength approach is crucial to identifying 
multiple-component systems. Indeed, the analysis of the 
radio morphology alone cannot distinguish between pairs 
of radio sources which are close in projection, or physi- 
cally connected radio components of the same source. In 
our sample there are 24 systems, whose radio morphol- 
ogy can be interpreted as multi-component radio sources. 
Their radio contours are plotted over the R image in Fig. 

El 

We perform the likelihood ratio analysis described in 
Section I3TT1 to look for a possible counterpart of each sin- 
gle radio component. For seven such systems (panels (d), 
( n )> (p)> ( s )> ( u ) in Fi &- 0), we find highly re- 
liable counterparts for each component. Therefore, we 
claim that they represent distinct sources. The other 17 
are confirmed to be multiple-component systems. They 
have extended radio emission, in most of the cases char- 
acterized by a core (not always visible in the radio) and 
radio lobes. The radio lobes have usually comparable 
radio power and are not associated with any optical or 
IR counterpart. There are some cases where we cannot 
exclude a possible contribution to the radio flux den- 
sity from a superimposed unrelated object like in sources 
with RID 38, 73, 209, 283, and 647. The complete list 
of these system with the properties of each radio compo- 
nent is given in N. Miller et al. (2012, in preparation). 
In Appendix 151 we discuss some peculiar sources. 

3.4. Revisited associations 

In section 13. 1\ we described our method to select the 
optical and IR counterpart for the radio sources using a 
likelihood ratio technique. Visually inspecting the results 
of the identification process, we confirm the association 
obtained following this procedure in 99% of the cases. 
In this section, we describe the reasons why we revisited 
the counterpart association for some peculiar sources. In 
these cases, the most likely counterpart has a reliabil- 
ity under our threshold and was therefore rejected. The 
two main reasons for this are: (i) the radio source is in 
a crowded field and therefore all the possible counter- 
parts have low reliability (RID 70, 407, 417, 458, 561, 
and 797). We based our choice of the counterpart on 
the radio morphology and on the overall object proper- 
ties in the various bands (as an example, see the notes 
on RID 407 in Appendix IB]) . (ii) The radio source is ex- 
tended, hence the exact position of the radio emission is 



not well determined or does not correspond to the opti- 
cal/IR peak emission (RID 407, 420, 521, 804, 828, and 
830). As a consequence, there is an offset between the 
position of the radio source and the counterpart that has 
therefore a low reliability. 

3.5. Estimation of spurious associations and 
comparison with the cross- correlation method 

For each auxiliary catalog, we estimate the rate of spu- 
rious associations by randomly shifting the position of 
the radio sources and computing again the reliability for 
all the possible counterparts. We apply only shifts be- 
tween 5 and 15 arcsec in order not to exceed the field 
coverage. We then compute the likelihood ratio value for 
each one of the shifted sources using equation [TJ where 
q(m) and n{m) are the probability distr ibutions derived 
for the original catalog (see Section I3TT]) . The same reli- 
ability threshold of 0.6 is adopted. The average fraction 
of false association over 50 different shifts is reported in 
Table [2] for each auxiliary catalog. We find spurious frac- 
tions from 3% up to 5% for the deep optical catalogs. In 
the case of H-GNS the fraction is very low but it could 
be underestimated due to the small area covered by this 
catalog and the consequent small statistics. 

We compare these results with the number of spuri- 
ous associations obtained using a simple cross-correlation 
method. The matching radius chosen for this test is 
equal to the radius which includes 90% of the counter- 
parts identified with the likelihood-ratio method. These 
radii are listed in column (4) of Table O We find that 
the two methods identify a similar fraction of sources for 
the optical catalogs and a somewhat lower one for the 
IR catalogs. We estimate the fraction of spurious asso- 
ciations similar to what has been done for the likelihood 
method, namely shifting the radio catalog with respect 
to the auxiliary one. We used the same set of displace- 
ments as in the previous case. As shown in Table [2j 
our likelihood ratio technique is generally less affected 
by spurious contamination especially when applied to the 
deep optical catalogs. In particular, in the cases of the 
U-VIMOS and v-GEMS catalogs the spurious fraction 
exceeds 10% with the cross-correlation method. If we 
decrease the searching radius from 0.7 to 0.5 arcsec for 
all the optical catalogs, the fraction of false counterparts 
becomes lower (8%, 7%, 3%, for the U-VIMOS, v-GEMS 
and R-WFI catalogs, respectively), but we also miss a 
significant fraction of real identification (19%, 23% and 
14%, respectively). For NIR and MIR catalogs, the two 
methods are almost equivalent. We obtain slightly higher 
fractions of fake associations with the likelihood ratio 
technique but with the cross-correlation method we miss 
a larger number of counterparts. This is mainly due to 
the lower source surface density with respect to the opti- 
cal catalogs. We note that the shift-and-rematch method 
tends to overestimate the number of false matches as it 
ignores the fact that there are a lar ge fraction of the 
sources that do have c ounterparts (see Broos et alJ 2007. 
l2011l : IXue et al.ll201lL for details.). Our estimates should 
therefore be considered as upper limits. Since this ef- 
fect is the same both for the likelihood and the cross- 
correlation method, it does not affect our conclusions. 
Finally, we assume that the fraction of spurious associa- 
tion in the final catalog is equal to the weighted average 
of the spurious fraction of each catalog, using the num- 
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Figure 5. Left panels: radio contours plotted over a 20x20 arcsec R band image. The RID is shown on the top right of each 
image. Other panels: cutouts of the images in the various bands centered at the position of the radio source (marked by a 
cross). Each image is 10 arcsec on a side. The corresponding catalog is indicated on the bottom left. The square indicates the 
position of the selected counterpart. The image from which the counterpart is selected has a black border. 



ber of counterparts selected from each catalog as weight. 
We then expect at most 4% spurious counterparts. 

3.6. Comparison with previous work 



The brighter sources of the 
already included in the radio 



Kellerm ann et al 



present catalog were 
catalog described in 



(120081). We h a ve co mpared the coun- 



terparts found in Ma inieri et~aT1 (|2008[ ) (M08, hereafter) 
with those selected in this work for the sources in com- 
mon. We find that in 90% of the cases the same counter- 
part is selected. For the remaining 10%, we select a dif- 
ferent counterpart compared to M08. In most cases the 
counterparts were identified in the optical band by M08 
while we find a more convincing counterpart in newly 
acquired MIR observations. For eight out of twelve ob- 
jects that were unidentified or under the threshold in 
M08, we now have a r eliable counterpart . One example 
is RID 625 (ID in the iKellermann et all (|2008f ) catalog 
KID=202) that was in an empty field in M08 while we 
identify it at 24 /im. The criteria to identify the best 
counterpart candidate of a radio source are slightly dif- 
ferent between the two works. In M08 a likelihood ratio 
threshold of 0.2 was adopted, while we use a cut in re- 
liability to take into account the presence of multiple 
counterpart candidates of the same radio source. Our 
reliability threshold of 0.6 is a bit more conservative (it 
correspond to a LR threshold of ~ 0.3) and it is aimed 



at reducing the number of spurious identifications. An- 
other difference between the two works is that in M08 
the best counterpart was chosen according to an a priori 
ranking of the auxiliary catalogs. The priority was set 
according to the depth of the optical/NIR survey and 
to the wavelength. Given the larger number of auxiliary 
catalogs used in this work we use a different approach: 
between catalogs in the same wavelength range (optical, 
NIR and MIR), we select the counterpart with the high- 
est reliability (see Section I3T]) . This allows us to fully 
exploit the information given by the probability distribu- 
tions obtained with the likelihood ratio technique. More- 
over, we minimize the number of tentative associations 
selected in a high priority optical/NIR catalog yet with 
a reliability just above the threshold. As a consequence, 
with our new approach 90% of the counterparts have a 
reliability greater than 0.96, in contrast to 0.83 obtained 
in M08 (see Fig. 4 in M08). Moreover, we observe a sig- 
nificant decrease in the average separation between the 
radio source and its counterpart compared to M08. In 
this work, we find 90% of the counterparts within 0.7 arc- 
sec around the radio position, which is about half of the 
radius found in M08 (see Fig. [7]). This is partially due to 
the change in r esolution of the two rad io surveys (from 
3.5" x 3.5" in IKellermann et all ([20081 , to 2.8" x 1.6" 
in this work). 

In summary, using the new MIR imaging in the E- 
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Figure 6. Cutouts of radio sources with complex radio morphology. The majority of them (17/24) are confirmed as multiple- 
component system, while for the remaining we find reliable counterparts for each radio component (see Section T3. 3 1 for details). 
Radio contours are plotted over the R band WFI image. The scale of the cutout is given in arcsec on the bottom right and the 
RID on the top left. 

CDFS area, we reach the same fraction of identification 
(95%) as in M08, although we adopt a higher reliability 
threshold and a larger fraction of the radio sources is in 
the outer part of the field, where the multi-wavelength 
coverage is poorer. There is only one source previously 
identified that we now consider unidentified: RID 101 
(see details in Appendix [Bj . 



4.1. 



4. REDSHIFT ASSOCIATIONS 
New VIMOS spectra and redshifts 



We acquired new optical spectra with the Visi- 
ble Multi-Object Spectrograph (VIMOS; iLe Fevre efaTI 



2003) at Very Large Telescope (VLT). We carried out 
one pointing in the central region of the VLA survey 
using the low-resolution (LR) blue grism (R=180, dis- 
persion=5.7 Apixel -1 ) that covers the wavelength range 
3700-6700 A. The total exposure time of five hours 
was set to identify faint optical counterparts to a lim- 
iting point-source magnitude of R~ 25. The mask was 
designed with the VIMOS M ask Preparation Software 
fVMMPS: lBottini et al1l2005l ) that optimizes the slit as- 
signments based on our input catalog. We observed a 
total of 32 VLA sources. The data were reduced using 
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Table 4 

List of the spectroscopic catalogs considered in this work. The instrument used to obtain the spectra is reported in column (2) and details 
on the observation and data reduction can be found in the references given in the first column. The label in column (3) are used in the 
final catalog to identify the source for the spectroscopic redshift. The last column reports the number of spectroscopic redshift adopted 

from each catalog in this work. 



(1) 


(2) 


(3) 


(4) 


Reference 


Instrument 


Label 


Number of z-spec adopted 


This work 


VIMOS 


P81 


13 


Szokolv et al. C2004) 


FORS1/FORS2 


S04 


38 


Vanzella et al. (2008) 


FORS2 


FORS 


20 


Silverman et al. (2010) 


VIMOS 


VJB 


37 


Silverman et al. (2010) 


VIMOS 


P80 


3 


Silverman et al. (2010) 


Keck 


K07 


32 


Silverman et al. (2010) 


Keck 


K08 


18 


Treister et al. (2009) 


VIMOS 


T09 


20 


Balestra et al. (2010) 


VIMOS-LR 


VLR 


24 


Balestra et al. (2010) 


VIMOS-MR 


VMR 


23 


Le Fevre et al. (2004) 


VIMOS 


VVDS 


19 


Ravikumar et al. (2007) 


VIMOS 


R07 


12 


Szokolv et al. (2004) 


FORS1/FORS2 


S04F 


1 


J. Kurk et al., 2012 (Submitted) 


FORS2 


GMASS 


4 


Norris et al. (2006) 


2dF 


N06 


10 




distance (arcsec) 



Figure 7. Cumulative distribution of the separation between 
a radio source and its counterpart for the radio catal o g pre - 
sented in this work (solid line) and in Maini eri et al.l (2008) 
(dot-dashed line). The vertical lines mark the separation 
within 90% of the counterparts are located. 

the VIM OS Interactive Pipelin e and Graphical Interface 
(VIPGI; iScodeggio etaLl I2005D . and the redshifts were 
estimated using the EZ0 software that cross-correlates 
each spectrum with a template spectrum, and via visual 
inspection to validate the result. We derived a spectro- 
scopic redshift for 13 VLA sources for which previously 
we had only a photometric redshift estimate. The spec- 
tra of these 13 radio sources are shown in Appendix [Al 

4.2. Spectroscopic redshifts 

9 http: / /cosmos . iasf-milano. inaf . it /pandora/EZ . html 




Redshift 

Figure 8. Top panel: fraction of photometric redshifts as a 
function of redshift. Bottom panel: total (empty) and spec- 
troscopic (filled) redshift distributions. The plot is cut at 
z = 4 for a better view of the range where both photo-z and 
spec-z are available. 

Many spectroscopic campaigns have been conducted 
in the E-CDFS. A complete reference list for those used 
in this work can be found in Table HI We combine the 
publicly available redshifts with our own newly acquired 
spectra (see Sec. 14. ip . We assign a quality flag (QF) to 
each redshift by mapping the ones in the original catalogs 
to a uniform scale. We use QF=3 to indicate a secure 
redshift, QF=2 for reasonable redshift, and QF=1 for 
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Figure 9. Radio power as a function of redshift. Circles 
indicate spectroscopic redshifts while crosses indicate photo- 
metric redshift s. The solid line shows the radio flux density 
limit of the VLA survey. 



tentative redshift or for single line detection. The co- 
ordinates of the counterparts, identified as described in 
Sec. 13. 1( are cross-correlated with the sources reported 
in the spectroscopic catalogs within 0.2 arcsec. Such a 
small radius is chosen to minimize the confusion with 
nearby sources. We find a spectroscopic redshift for 274 
sources. If a source has a match in more than one spec- 
troscopic catalog, we verify the consistency between the 
corresponding redshifts. For 22 sources, multiple spec- 
troscopic redshifts differ by more than 0.1. In all but 
three of these cases, the QF of the spectroscopic mea- 
surements allows us to select the highest quality z. For 
sources RID= 83, 569, and 706, where the QF in the 
various spectroscopic catalogs are equivalent, we visually 
checked the spectra to select the more reliable redshift 
value. The spectroscopic redshift associated with each 
radio source is reported in the column 7 of Table [5j with 
the QF, and the reference in columns 8 and 9, respec- 
tively. In summary, 33% of the radio sources with coun- 
terpart have a spectroscopic redshift, 74% of which are 
secure redshift (QF=3), 18% have QF=2 and 8% have a 
tentative redshift measurement (QF=1). 

4.3. Photometric redshifts 

In order to increase the redshift completeness of our 
sample, we also use photometric redshift estimates. 
We use the pho tometr ic redshift cata l og co mpiled by 
ILuo et all (|2Q1QD and iRafferty et all (|2011[ ). These 
redshifts are based on a large number of photomet- 
ric b a nds: the COMBO -17 optical catalog ([Wolf et al.l 
l200l iWolf et all l2008h. the GOODS-S MUSIC cata- 
log (Grazian et al. 2006), the MUSYC BVR-detected 
catalog (Ga wiser et a l. 2006), the deep G OODS-S VI- 
MOS U-band catalog (jNonino et al.ll2009[ ). the GALEX 
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Figure 10. Normalized redshift distribution for the radio 
sources with an optical counterpart (top), with a NIR coun- 
terpart (middle) and with a MIR counterpart (bottom) . The 
mean redshift is increasing from 1 for optical identified sources 
to 2.5 for the MIR ones (see Section O]). 

Data Rele ase 43, the MUSYC near-infrared catalogue 
dTavlor et aflj2QQ9|) : and the SIMPLE mid-infrared one 



(Da men et al 



l201l( ). Starting from this photometric 
data set, the publicly available Zuric h Extragalactic 
Baye sian Redshift Analyzer (ZEBRA; iFeldmann et al.l 
2006) code was used to derive photometric redshifts via 
a maximum likelihood technique. The set of templates 
used includes: 259 galaxy templates constructed from 
PEGASE stellar population synthesis models, a set of 
hybrid (galaxy +AGN) templates an d ten empirical AGN 
te mplates fromlPolletta et all (120071). We refe r the reader 
to ILuo et all (j2010f ) and IRafferty et al.l (|2011f ) for a more 
detailed description of the procedure adopted to estimate 
photometric redshifts. We cross-correlate the photomet- 
ric redshift catalog with the radio source counterparts, 
selected as described in 13.21 Given the high background 
surface density distribution of the photometric catalog, 
we adopt a matching radius of 0.2 arcsec. This way we 
minimize the risk of associating to the radio counterpart 
the redshift of a nearby source. We find 623 matches out 
of the 839 identified radio sources. The mean separation 

10 http://galex.stsci.edu/GR4/ 
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between the radio counterparts and the corresponding 
photometric redshift coordinates is 0.03 arcsec. For the 
remaining 216 objects we consider three other compila- 
tions of photometric redshifts : the MUSYC-E-CDFS cat- 
alog (Cardamone et al. 2010), the GOODS-MUSIC cat- 
alog (Sa ntini et aTll2QQ9ft. and the K-selected MUSYC 
catalog (Tayl or et al.l l2QQ9h . The latter is based on a 
NIR selected catalog, and therefore is particularly useful 
to assign a redshift to radio sources with no counter- 
part in the optical. We find 13, 4, and 35 additional 
redshifts, respectively. In Summary, we associate a pho- 
tometric redshift to 673 (80%) out of the 839 identified 
radio sources. 

For the sub-sample with available secure spectroscopic 
redshift (QF=3), we compute the normalized median ab- 
solute scatter, 

. , ,. f \Az — mediant Az)\\ . N 
o~nmad = 1-48 x median I I , (5) 

\ 1 H~~ %spec J 

where Az = z p hot — z speci which is an estimate of the 
quality of the photometric redshift whic h is less sensitive 
to ou tliers than the standard deviation ([Brammer et al.l 
2008). We find o~nmad = 0.01 that is comparable, 
and even slightly better, to what found for the same in- 
dicator in the photometric redshift catalogs considered 
(ISantini et al . 2009; Taylo r et al1l2QQ9b ICardamone et all 
l201Qt [Raffertv et al.l 1201 J ). Therefore, we conclude that 
the accuracy of the photometric redshifts for our radio 
selected sample is comparable to what estimated for the 
overall population in the photometric catalogs. We note 
that this indicator assumes that the spectroscopic sub- 
sample is representative of the full sample. This assump- 
tion is likely not entirely true and consequently it gives 
an overestimat ion of the accura cy of the photo-z for the 
whole sample. iLuo et aTl (j201Q[ ) estimated that the un- 
certainties on the photometric redshifts for the full sam- 
ple are a factor of six higher. We also note that the 
photomet ric redshift errors giv en in column 6 of Table [5] 
from the Ra ffertv et al.l (120111) catalog are known to be 
underestimated (see ILuo ^Tall (jl010)). A more realis- 
tic estimate is given by the ctnmad parameter, which is 
around 6% as discussed above. 

4.4. Redshift distribution 

In the case where both spectroscopic and photometric 
redshifts are available, we use the spectroscopic one if QF 
> 2 and the photometric redshift otherwise. For spec- 
troscopic redshifts with QF=1 we find ctnmad — 0.28, 
which we interpret as an indication of the poor quality of 
these spectroscopic redshifts. Combining spectroscopic 
and photometric information, we assign a redshift to 678 
objects, 81% of the radio sources with counterpart (252 
spectroscopic redshifts and 426 photometric redshifts). 
This fraction underestimates the redshift completeness 
of our sample since in the outermost part of the field 
there are no redshift measurements available. There- 
fore, we restrict our redshift distribution analysis to the 
sources in the area covered by the photometric redshift 
catalog^]. This region is plotted with a dot-dashed line 
in the right panel of Fig. [2] The number of radio sources 

11 Only three sources outside this region have a spectroscopic 
redshift. 



included is 779, and 87% of them have a redshift. The 
total redshift distribution is plotted in Fig. [8j where 
the filled histogram represents the distribution of sources 
with spectroscopic redshifts. The top panel shows the 
fraction of photometric redshifts. We note that photo- 
metric redshift measurements become increasingly im- 
portant at higher redshifts where optical spectroscopic 
observations become more challenging. The mean red- 
shift for the whole sample is 1.1 and the median is 0.9. 
Fig. [9] shows the radio power as a function of redshift 
for sources with either spectroscopic (circle) or photo- 
metric (cross) redshift. The solid line represents the 
flux density limit of the survey. If we divide the radio 
sources based on their identification band, we observe 
an increase in the mean (median) redshift from 1.0 (0.8) 
for the optical identified sources to 2.5 (2.1) for the MIR 
ones (see Tableland Fig. [10]). The statistical signifi- 
cance of the different redshift distributions is examined 
with the Kolmogorov-Smirnov (K-S) test. The difference 
between the optical and NIR distributions and that be- 
tween the optical and MIR distributions are confirmed 
with a significance level > 99%. For the NIR and MIR 
redshift distributions the K-S test gives a significance of 
99%. Both spectroscopic and photometric redshift esti- 
mates become more challenging moving to high redshift 
objects and this is the reason why the fraction of sources 
with a redshift estimate drops from 88% for the optical 
identified sources to 25% for the MIR ones (see Table [6]). 

5. X-RAY COUNTERPARTS 

Chandra has imaged in the X-ray the area of the E- 
CDFS as part of two different programs. The first is 
a 250 ks exposure obs ervation that cover s almost the 
whole field (0.28 deg 2 ) (jLehmer et al.ll2005D . The survey 
reaches sensitivity limits of 1.1 x 10 -16 and 6.7 x 10 -16 
ergs cm -2 s _1 for the soft (0.5 - 2.0 keV) and hard (2 
- 10 keV) bands, respectively. The second set is a much 
deeper 4 Ms Chandra observation covering only the cen- 
tral part of the field (~ 0.1 deg 2 ). The on- axis flux limits 
are 9.1 x 10 -18 for the soft band and 5.5 x 10 -17 ergs 
cm -2 s _1 for the hard band (|Xue et al.ll2011f ). We cross- 
correlated the radio source catalog with the X-ray ones. 
Due to the low surface density of X-ray sources a simple 
positional match is almost equivalent to the likelihood 
ratio technique. The searching radius was set to three 
times the sum in quadrature of the errors on the radio 
and X-ray positions. In case of multiple counterpart can- 
didates, we selected the one closest to the radio source 
position. We find 129 radio sources with X-ray detection 
from the 4 Ms Chandra catalog, and 99 sources from 
the 250 ks catalog. Combining the two lists, we have 
X-ray detection for 25% of our radio sources. Their flux 
in the soft and hard band are reported in columns 10 
and 1 1 of Table HI We re fer the reader to Section 4.2 of 
IVattakunnel et aTl (j2012[ ) for a description of the proper- 
ties of the radio sources with X-ray counterpart. In the 
following section we focus our attention on the sources 
for which we obtain only upper limits on their X-ray flux. 

5.1. Average X-ray properties of radio- only detected 
sources 

The majority (75%) of our radio sources has no X-ray 
counterpart. Even in the region covered by the 4 Ms 
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Table 5 

Main characteristics and redshift information of the radio source counterparts. A description of the table content is given in Section [6] 
(This table is available in its entirety in a machine-readable form in the online journal. A portion in shown here for guidance regarding its 

form and content.) 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) (8) 


(9) 


(10) 


(11) 


(13) 


RID 


R mag 


Ks mag 


3.6/21X1 mag 


best-z 


photo-z 


spec-z QF 


reference 


S 0.5-2fceV 


S 2-10fceV 


XID 




(AB) 


(AB) 


(AB) 










(erg cm - ^ s_ "*") 


(erg cm - ^ s_ "*") 





361 


25.92 ± 


0.22 


22.02 


± 


0.18 


21 


02 


± 0.02 


1.39 


1 


39 


362 


20.40 ± 


0.00 


18.47 


± 


0.01 


18 


46 


± 0.00 


0.61 





60 


363 












21 


33 


± 0.30 








364 


25.57 ± 


0.14 








22 


50 


± 0.07 


3.02 


3 


02 


365 


22.78 ± 


0.02 


19.60 


± 


0.03 


19 


29 


± 0.00 


0.73 





73 


366 


22.91 ± 


0.02 


19.98 


± 


0.04 


19 


73 


± 0.00 


0.64 





64 


367 


21.17 ± 


0.00 


19.16 


± 


0.02 


19 


16 


± 0.00 


0.60 
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Table 6 

Redshift distribution of the counterparts divided according to their wavelength selection. The number in column (1) is the number of 
radio sources whose counterpart is found in the group of catalogs of the corresponding row. 





(1) 


(2) 


(3) 


(4) 


(5) 




Identified sources 


With z 


Fraction 


Mean z 


Median z 


Optical 


652 


575 


88% 


1.0 


0.8 


NIR 


76 


75 


98% 


1.6 


1.5 


MIR 


111 


28 


25% 


2.5 


2.1 



Chandra obser vation, the fraction of identified sources 
is only - 60% (jVattakunnel et al.ll2012l ). For the radio- 
only sources we perform aperture photometry on the X- 
ray images at the position of the radio source. The X-ray 
detected sources are masked and replaced with a Poisso- 
nian background based on the value of the measured local 
background. The photometry is done separately in the 
soft (0.5-2 keV) and in the hard (2-10 keV) bands. To 
derive the average properties of these objects, we perform 
a stacking analysis of the Chandra images. In particular, 
we stack separately sources with counterparts selected 
from an optical catalog, from a NIR catalog and from a 
MIR one. The net counts obtained in both hard and soft 
X-ray bands are reported in Table The detection for 
the optical selected sources is highly significant in both 
bands, while for the MIR ones it is only marginal in the 
soft band. The NIR selected sources have marginal de- 
tection only in the soft band and are not detected in 
the hard band. For each group, we evaluate the aver- 
age hardness ratio defined as HR = (H — S)/(H + 5), 
where H and S are the total net counts in the hard and 
soft band, respectively (column 4 in Table [7]). In par- 
ticular we note that MIR selected sources have a hard 
hardness ratio, HR = 0.4 ± 0.2, supporting the hypoth- 
esis that these objects are obscured sources. This HR 
value corresponds to an effective X-ray photon indices 
r = 0.07t°;H 

We also split the sample of X-ray undetected sources in 
radio power bins to investigate if there are any changes in 
the average X-ray spectral properties as a function of ra- 
dio power. We consider only sources with z < 1.5, where 



we have a more uniform distribution in radio power, from 
10 20 to 10 27 W Hz -1 . The radio power bins, the net 
counts, and the HR are reported in Table [71 We find a 
roughly constant value of HR and therefore no signifi- 
cant change in the average X-ray spectral properties as 
a function of radio power. 

6. CATALOG DESCRIPTION 

In this section we describe the catalog containing the 
results of the optical and IR counterpart identification 
process. The information is divided into two tables. In 
Table [3] we include the radio data from N. Miller et al. 
(2012, in preparation) that were used in this work and 
the results of the identification process. In Table [5] we list 
the main characteristics of the optical or IR counterpart, 
the redshift information and the X-ray data. The catalog 
columns are organized as follows. In Table [3j 

• (1) Identification number of the radio source (RID). 

• (2) and (3) Right ascension and declination of the 
radio source. 

• (4) Radio flux density and la error in /iJy. 

• (5) Signal to noise ratio. 

• (6) and (7) Right ascension and declination of the 
counterpart. 

• (8) Reliability of the association^. 

12 Sources whose identification has been revisited (see Section 
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Table 7 

Net counts and HR for radio sources without X-ray detection. The number in column (1) is the number of radio sources which have a 

optical/IR counterpart but not detected in the X-ray. 





(1) 


(2) 


(3) 


(4) 




Number of sources 


0.5-2 keV cts 


2-10 keV cts 


HR 


Optical 


426 


676 ± 68 


374 ± 97 


-0.3 ± 0.1 


NIR 


59 


78 ± 24 






MIR 


86 


79 ± 33 


227 ± 50 


0.4 db 0.2 


logP r < 22.9 


122 


231 ± 36 


162 db 51 


-0.2 ± 0.2 


22.9< logP r <23.5 


121 


194 ± 35 


98 ± 50 


-0.3 ± 0.2 


logP r > 23.5 


97 


120 ± 31 


52 ± 43 


-0.4 ± 0.4 



(9) Distance between the radio source and the 
counterpart in arcsec. 

(10) Catalog from which the counterpart is selected 
(see Table [Q). 

Table El 

(1) Identification number of the radio source (RID). 

(2) R-band AB magnitude of the counterpart from 
the WFI catalog and associated error. 

(3) K-band AB magnitude of the counterpart from 
the MUSYC catalog and associated error. 

(4) Flux density at 3.6/im of the counterpart from 
the SIMPLE catalog and associated error. 

(5) Best redshift of the counterpart: spectroscopic 
if QF > 2, photometric otherwise. 

(6) Photometric redshift with upper and lower 68% 
confidence level. 

(7) Spectroscopic redshift. 

(8) Quality flag (QF): 3 for secure redshift, 2 for 
resonable redshift and 1 for one line detection or 
tentative redshift. 

(9) Source of the spectroscopic z (see Table HJ. 

(10) X-ray soft band flux (0.5-2.0 keV) and associ- 
ated error. 

(11) X-ray hard band flux (2-10 keV) and associ- 
ated error. 



(from IXue et al 



(|2005h 



(12) X-ray ID 
<1000, from iLehmer et al.l 
[ID-Lehmer+1000]) 

7. DISCUSSION AND CONCLUSIONS 



(|2oTTh if 

if >1000 



We have presented the optical and IR counterparts of 
the radio sources in N. Miller et al. (2012, in prepara- 
tion) catalog. The results are collected in a new catalog0 
containing the counterpart data and the redshift infor- 
mation. A detailed characterization of the physical prop- 
erties of these sources will be presented in M. Bonzini et 

13. 4|) or for which the counter part position has been extracted from 
the IRAC image (see Section [3?2]) have Rel = 9. 

13 The catalog is available in ASCII format in the on-line mate- 
rial. 



al. (2012, in preparation). This work has demonstrated 
the difficulties in, and the requirements for, the identifi- 
cation of the sub-mJy radio population. The main results 
of our analyses are as following: 

1. Importance of multi-wavelength observations. We 
identify the counterparts for a high fraction (95%) 
of radio sources. In order to reach such a com- 
pleteness it is necessary to include not only op- 
tical observations, but also near and far infrared 
data. Optical surveys alone, even in the deepest 
fields, allow us to identify only ~ 70% of the radio 
sources. With just MIR observations the fraction 
rises to 86%, but it is by only combining the in- 
formation from all wavelengths that we reach 95% 
completeness. The mult i- wavelength coverage is 
also important to obtain a high redshift complete- 
ness. Indeed, only 31% of our radio sources have 
spectroscopic information, while the majority have 
a photometric redshift. 

2. Importance of the counterpart analysis to confirm 
multiple- component radio sources. In this work we 
have found many examples that show the impor- 
tance of the combination of radio and optical/IR 
data to correctly identify multiple-component ra- 
dio system. In many cases, sources whose radio 
morphology suggested a complex radio structure 
(e.g., KID 114) have been identified as indepen- 
dent sources. The opposite case is represented by 
source RID 73. Here, we conclude that the radio 
emission is associated with a single compact radio- 
lobe source with a single optical counterpart. 

3. Comparison between likelihood rat io and cross- 
correlation methods. In Sec 13. 5[ we compare 
the likelihood ratio technique with the positional 
matching method. This work has shown that the 
latter is hardly applicable to deep optical surveys 
since it leads to a large fraction of spurious matches 
(~ 10%). With our technique instead the rate of 
spurious matches is lower due to the exploitation of 
the information given by the probability distribu- 
tion of background sources in the optical catalogs. 
We have also shown that to reach the same level of 
spurious contamination with the cross-correlation 
method the fraction of identified sources decreases 
by ~ 18%. At longer wavelengths, i.e. in the NIR 
and in the MIR, the differences of the two methods 
are negligible. We find a comparable fraction of 
expected spurious counterparts and a similar com- 
pleteness. This is mainly due to the lower back- 
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ground surface densit y of objects in the auxiliary 
catalogs. In section 13.41 we point out two cases 
where the likelihood ratio method can fail in iden- 
tifying the correct optical counterpart, that is the 
presence of many close sources around the radio 
object and extended radio emission. We note that 
with the cross-correlation method these problems 
are even more severe; in crowded regions it selects 
the closest source, regardless of its properties. In 
case of extended radio sources, there could be an 
offset between the peak of the emission at different 
wavelengths larger than the searching radius. In 
this case, the cross-correlation method would not 
be able to give any counterpart candidate. 

4. Comparison with MO 8 work. Compared with the 
sample studied in M08, our sample is about 3 times 
larger and most of the new sources have a low ra- 
dio flux density or lie at the edges of the E-CDFS. 
That makes their identification more challenging. 
However, more and deeper catalogs are now avail- 
able in the E-CDFS and, using these data, we were 
able to reach the same identification completeness 
as in M08 for the new sample. We find general 
agreement between the counterparts found in the 
two works for the radio sources in common. The 
main improvement is a more reliable identification 
in particular of the optically faint radio sources, 
obtained by adopting a stricter acceptance crite- 
ria and giving more importance to the IR selected 
catalogs (Section 13. 6p . 

5. Importance of MIR observations to find the coun- 
terpart of high redshift or heavily obscured radio 
objects. Some radio sources (12%) have a reli- 
able counterpart only in the catalogs based on the 
Spitzer data. These sources are particularly in- 
teresting since they are the best candidate high- 
z objects. In Sec. 14. 4| we describe the redshift 
distribution of the radio sources divided accord- 
ing to their identification band. Indeed, we find 
a clear trend for sources identified at longer wave- 
lengths to have higher redshifts, as shown in Fig. 
flQl Moreover, the stacking analysis of the X-ray 
images of the MIR selected sources, has revealed 
that they tend, on average, to have hard X-ray 
spectra (HR = 0.4). This supports the idea that 
they are obscured sources. 
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APPENDIX 
NEW VIMOS/VLT SPECTRA. 

Newly acquired VIMOS spectra for the counterparts of 13 VLA sources (see Sec. 14. ip . We used the low-resolution 
blue grism (R=180), with a total exposure time on source of 5 hours. Each plot of Fig. [TT1 shows the spectra and the 
position of the main spectral features. We make bold the names of the lines actually used to identify the redshift. The 
labels report the source RID, the redshift, and the corresponding quality flag (QF). 




400 450 500 550 600 650 400 450 500 550 600 650 

Observed wavelength [nm] Observed wavelength [nm] 



o 




400 450 500 550 600 650 400 450 500 550 600 650 

Observed wavelength [nm] Observed wavelength [nm] 




400 450 500 550 600 650 400 450 500 550 600 650 

Observed wavelength [nm] Observed wavelength [nm] 



Figure 11. New VIMOS/VLT spectra of the radio sources. 
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Figure 11. (Continued.) 
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NOTES ON INDIVIDUAL SOURCES. 

RID 73 (KID 14): there are two possible counterparts for the radio lobes but we believe that they are both 
associated with the bright central galaxy. One strong indication for this is the strength of the radio flux density, 
which at 40 mJy is reasonable for a compact double lobed radio galaxy. If separate sources, they would have to 
both be strong radio AGN very close in projection on the plane of the sky. 

RID 80-85 (KID 18B-18A): they were considered as radio lobes in Kelle rmann et all (|2008[ ). but they appear to 
have two different counterparts. Moreover, they have very different flux densities which supports the idea that 
they are not related to the same source. Finally, there is no good candidate for a single radio core. 

RID 101 (KID 23): double lobed source whose core was tentatively identified in M08 (KID 23) with a faint 
(R-band magnitude ~ 26) galaxy at z=0.999 from the COMB017 catalog. We think that this association is 
very unlikely especially because this possible counterpart is not detected at any longer wavelength. From the 
radio contours we believe that it is a classical radio gala xy. Therefore, we exp ect for such a galaxy a correlation 
between the K-band magnitude and the redshift (e.g.. iLilly fc Longair| [l984). Since we do not detect it in the 
K-band we think this source may be a high redshift object. This hypothesis is also supported by the presence of 
a possible counterpart at 24 micron in the FIDEL catalog. Unfortunately no photometric redshift is available. 

RID 209 (KID 48): V-shape radio source at redshift 1.3. Given the radio morphology, we considered the 
possibility that this is a head-tail radio source due to the interaction with the intracluster medium (ICM) in a 
high-z galaxy cluster. Any cluster, or large group, with ICM density sufficient to bend the radio jets, would have 
been clearly detected in the X-ray too. However, we do not observe it in the 4 Ms Chandra image. Also the 
redshift distribution of the sources in the region around RID 209 does not show any hint of clustering. Therefore, 
we think is unlikely that the V-shape of this radio source is due to the interaction with the ICM. There is instead 
a possible contamination to the flux density of one of the two radio lobes (A) from a superimposed galaxy. 

RID 283 (KID 73): double- lobe source. There are three optical sources in the region of the radio emission but 
all have reliability under the threshold (< 0.6) and we believe that none of them are associated with the radio 
source. The counterpart of the core is identified with an object detected in the Ks-MUSYC catalog. 

RID 308 (KID 80): bright radio source with possibly one or two lobes. However, the quality of the radio image 
in this region is not good. 

RID 360 (KID 97): powerful single component radio source that was not identified in M08. It has a radio flux 
density of 1.38 mJy. Although we use deeper catalogs, still we are not able to identify it in any band. The 
cutouts of this source are shown in Fig. [5] and they are all blank field. The 5a detection limit for each band is 
given in Table [TJ Moreover, this source is in the region covered by the 4Ms Chandra observation but it has no 
X-ray detection. We can therefore put an upper limit on its X-ray flux of 9.1 x 10 -18 and 5.5 x 10 -17 ergs cm -2 
s _1 for the soft and hard band, respectively. 

RID 403-406-410-412 (KID 114): this group of sources was at first interpreted as a tailed radio source (see radio 
contours). But since we find a clear counterpart for three of them, we consider these sources as independent. 
Source 410 is unidentified. 

RID 407 (KID 113): bright and extended double lobed source. Close to the core position there are many optical 
sources. In particular there is a 21 K-band magnitude galaxy 0.5 arcsec away from the expected core position 
that was automatically selected by our method. We consider t his association spurious since it would impl y 
that this source is far from the K-z relation for radio galaxies (e.g lLillv fc Lo ngair 1984; De Breuck et al.l f2002). 
Therefore, we manually corrected the identification by associating this radio source to a bright elliptical galaxy 
(K-band mag=18.4) 2 arcsec away from the centroid of the radio image. 

RID 500 (KID 148): this is a complex radio source. We identify a clear core and two radio lobes (KID 148 A and 
148B). There are two other components, KID 148D and 148E, possibly associated with this source. We found a 
secure counterpart for 148D and so we listed it as a separated source (RID 504). For the 148E component the 
only counterpart candidate is a faint galaxy at a =03:32:32.59, 5 =-28:03:15.4 with R-mag=23.7, but it is under 
our reliability threshold and therefore it probably remains unidentified. 

RID 848 (KID 260): the radio source is split into two components both associated with the same spiral galaxy. 
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